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Abstract. We prove that the critical points of the 3d nonlinear elasticity functional on shells of 
small thickness h and around the mid-surface S of arbitrary geometry, converge as h — *• to the 
critical points of the von Karman functional on S, recently derived in [10] . This result extends 
the statement in [16], derived for the case of plates when S C I 2 . We further prove the same 
convergence result for the weak solutions to the static equilibrium equations (formally the Euler- 
Lagrange equations associated to the elasticity functional). The convergences hold provided the 
elastic energy of the 3d deformations scale like h 4 and the external body forces scale like h 3 . 
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1. Introduction and statement of the main results 

Since the beginning of research in nonlinear elasticity, a major topic has been the derivation of 
lower dimensional theories, appropriately approximating the three dimensional theory on struc- 
tures which are thin in one or more directions (such as beams, rods, plates or shells). Recently, 
the application of variational methods, notably the T-convergence [5], lead to many significant 
and rigorous results in this setting [91 [7] . Roughly speaking, a T-limit approach guarantees the 
convergence of minimizers of a sequence of functionals, to the minimizers of the limit. However, it 
does not usually imply convergence of the possibly non-minimizing critical points (the equilibria) 
and hence other tools must be applied to study this problem. 

In this note, following works |13| ,I141[T6] in which beams, rods and plates were analyzed, we study 
critical points of the 3d nonlinear elasticity functional on a thin shell of arbitrary geometry, in 
the von Karman scaling regime. A T-convergence result in this framework was recently derived in 
[10J, providing the natural from the minimization point of view generalization of the von Karman 
functional [7] to shells. In analogy with the analysis done in [16] for plates, we now proceed to 
establish convergence of weak solutions to the (formal) Euler-Lagrange equations (jl.lip . as well 
as convergence of critical points of the 3d energy functionals ([l.ip . to the critical points of the 
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functional obtained in |10j . As pointed out in [1] Problem 5, in general it is still unknown whether 
these two definitions of equilibria are equivalent. 

We now introduce the basic framework for our results. We consider a 2-dimensional surface 
5 embedded in R 3 , which is compact, connected, oriented, of class C 1 ' , and with boundary dS 
being the union of finitely many (possibly none) Lipschitz curves. A family {S h }h>o of shells of 
small thickness h around 5 is given through: 

S h = {z = x + tn(x); x G 5, -h/2 < t < h/2}, < h< h . 

By n(x) we denote the unit normal to 5, by T X S the tangent space, and II(x) = Vn(x) is the 
shape operator on 5 (the negative second fundamental form). The projection onto 5 along ft is 
denoted by it. We assume that h < ho, with ho sufficiently small to have ir defined on each S h . 
To a deformation u G W 1,2 (S h , M 3 ) we associate its elastic energy (scaled per unit thickness): 

(1.1) E h (u) = \ I WNu). 

h Js h 

The stored energy density W : M 3x3 — > [0, oo] is assumed to be C 2 in a neighborhood of 50(3), 
and to satisfy the following normalization, frame indifference and nondegeneracy conditions: 

VFeM 3x3 Vi?£50(3) W(R) = 0, W{RF) = W(F), 

{L2) W(F)> Odist 2 (F, 50(3)) 

(with a uniform constant C > 0). Our objective is to describe the limiting behavior, as h — > 0, of 
critical points u h to the following total energy functionals: 

(1.3) J\u) = E\u) - i / f h u, 

h Js h 

subject to external forces f h , where we assume that: 

f h (x + tn) = hVe^f(x)det(ld + tU)- 1 , f £ L 2 (S, M. 3 ) and f f = 0. 

Above, e h is a given sequence of positive numbers obeying a prescribed scaling law. It can be 
shown [10] that if f h scale like h a , then the minimizers u h of fjl .3j) satisfy E h (u h ) ~ hP with 
/3 = aifO<o<2 and (3 = 2a — 2 if a > 2. Throughout this paper we shall assume that (3 > 4, 
or more generally: 

(1.4) lim e h /h A = k < +oo, 

which for ScM 2 corresponds to the von Karman and the purely linear theories of plates, derived 
rigorously in [?]• 

In our recent paper [10], the T-limit of l/e h J h has been identified in the scaling range cor- 
responding to (|l,4p . and for arbitrary surfaces 5. It turns out that the elastic energy scaling 
E h {u h ) < Ce h implies that on 5 the deformations vfr s must be close to some rigid motion Qx + c, 

and that the first order term in the expansion of Q T (u^ s — c) — id with respect to h, is an ele- 
ment V of the class V of infinitesimal isometries on 5 [T7] . The space V consists of vector fields 
V S W 2 ' 2 (S,R 3 ) for whom there exists a matrix field A € W^ 2 {S, M 3x3 ) so that: 

(1.5) d T V(x) = A(x)t and A(x) T = -A(x) Va.e. x e 5 Vr G T X S. 

Equivalently, the change of metric on 5 induced by the deformation id + hV is at most of order 
h 2 , for each V £ V. 
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When in (jl.4[) k = 0, the limiting total energy is given by: 

(1.6) J(V, Q) = ^J^Q 2 (x, (y{An) - AU) tan ) dx - j f • QV dx, VV G V, Q G 50(3). 

The first term above measures the first order change in the second fundamental form II of S, 
produced by V. The quadratic forms Q 2 (x, •) are given as follows: 

Q 2 (x, F tan ) = mm{Q 3 (F); (F - F) tan = 0}, Q 3 (F) = D 2 W(ld)(F, F). 

The form Q 3 is defined for all F G M 3x3 , while Q 2 {x, •) for a given x G S, is defined on tangential 
minors i^n °f such matrices. Both forms depend only on the symmetric parts of their arguments 
and are positive definite on the space of symmetric matrices [6j. In the weak formulation of the 
Euler-Lagrange equations of (|1.6h one naturally encounters the linear operators £3 and C 2 (x, •), 
defined on matrix spaces M 3x3 and R 2x2 respectively, given by: 

VFGM 3x3 Q 3 (F) = C 3 F:F and Q 2 (x, F tan ) = C 2 (x, F tan ) : F tan . 

For k > 0, the T-limit (which is the generalization of the von Karman functional [7] to shells), 
contains also a stretching term, measuring the total second order change in the metric of S: 
(1.7) 

J vK (V, B tan ,Q) = lJ s Q2 (x, B tan - ^(A 2 ) tan ) +^J s Q2 (s, (y(An) - AU) tan ) - j f ■ QV. 
It involves a symmetric matrix field Btan belonging to the finite strain space: 

B = cl^jsymVu/ 1 ; w h G W 1,2 (S, R 3 ) j. 

The two terms in (|1.7p correspond, in appearing order, to the stretching and bending energies 
of a sequence of deformations v h = id + hV + h 2 w h of S, which is induced by a first order 
displacement V G V and second order displacements w h satisfying lim^o symVw' 1 = Bt an ■ The 
crucial property of (|1.7|) is the one-to-one correspondence between the minimizing sequences u h 
of the total energies J h (u h ), and their approximations (modulo rigid motions Qx + c) given by v 
as above with (V, Bt an , Q) minimizing J vK , or (V, Q) minimizing J when k = 0. 

The purpose of this paper is to show that under the following extra assumption of [16] : 

(1.8) VFGM 3x3 \DW{F)\ <C{\F\ + l). 

also the equilibria (possibly non-minimizing) of (jl.ip converge to the equilibria of (jl.7p or (jl.6p . 
The definition of an equilibrium of the 3d energy J h may be understood in two different manners, 
corresponding to passing with the scaling e of a variation <p to outside or inside the integral sign. 
Namely, for a fixed h > 0, we may require that: 

(1.9) G W^OS* R 3 ) lim 1 (j h (u h + e<j) h ) - J h {u h )^j = 0. 
or that: 

(1.10) V/ G W l ' 2 (S h ,M?) f DW(Vu h ) :V(f> h = [ f h <f> h . 

Js h Js h 

The last condition is obtained by formal passing to the limit e — > under the integral sign in (jl.9p . 
Integrating by parts we also see that (|1.10p is the weak formulation of the following fundamental 
balance law pQ: 

(1.11) div [DW(Vu h )] + f h = in S h , DW{Vu h )n = on dS h , 

where the operator div above is understood as acting on rows of the matrix field DW(Vu h ). 
Whether the two definitions of equilibria fjl .91) and (ll.lOp are equivalent, even for local minimizers 
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(without assuming extra regularity, e.g. their Lipschitz continuity) is an open problem of nonlinear 
elasticity, listed by Ball as Problem 5 in pQ. 

It turns out that the main convergence result described below follows with either (|1.9p or (jl.lOp . 
The reason is that the difference between these two definitions (after an appropriate scaling), 
converges to with h, along particular sequences of variations (f> , which are however exactly the 
3d variations recovered from the the variations of the 2d functional J vK or J. 

Theorem 1.1. Assume \1.2) and (fi.gj) . Let u h G W 1,2 (S h , M 3 ) be a sequence of deformations, 
satisfying: 

(a) the equilibrium equations hi. 10) hold, 

(b) E h {u h ) < Ce h , where e h is the scaling with ( fi.^| ). 

Then there exist a sequence Q h G 50(3), converging (up to a subsequence) to some Q G 50(3), 
and c h G M 3 such that for the normalized rescaled deformations: 

y h (x + tn) = (Q h ) T u h {x + h/h tn) - c h 

defined on the common domain S h °, we have: 

(i) y h converge in W 1,2 (S h °) to it. 

(ii) The scaled average displacements: 

h fV2 

(1.12) V h (x) = ^f y h (x + tn) 



x dt 

-h /2 ' 



converge (up to a subsequence) in W 1,2 (S) to some V G V. 

(iii) h/Ve^ sym X7V h converge (up to a subsequence) in L 2 (S) to some Bt an G B. 

(iv) The triple (V, Bt an ,Q) satisfies the Euler- Lagrange equations of the functional J vK . That 
is, for all V G V with A = \7V given as in the formula ( f and all Btan & B> there 
holds: 



■ B t an — 0, 



)tan 



11.14) 



(1.13) f C 2 (x,B t an-^(A 2 )tan) 

J s 

- K^£ 2 (x,B tan - ^(A 2 ) tan ) ■ (AA) t 

+ ^ j Ci (x, (V(An) - AU)tan) : (V(In) - AIi) tan = J f-QV, 

When k = then the couple (V, Q) satisfies (7777]) f or a ^ V S V, which is the Euler- 
Lagrange equations of the functional |77j|). 

Theorem 1.2. Theorem M . 1\ remains true if in the assumption (a) the formal equilibrium equation 
&1.1(J\) are replaced by the critical point condition lil.fA ). 

We prove Theorem 11.11 in section [2] and Theorem 11.21 in section [3j In section [5] we derive the 
third Euler-Lagrange equation (after the first two (| 1 . 13[) and (|1.14p ). corresponding to variation in 
Q G 50(3). We first notice that the limiting Q necessarily satisfies the constraint of the average 
torque r{Q) = J s f X Qx dx being 0. The main difficulty arises now from the fact that the 
variations must be taken inside 50(3) in a way that this constraint remains satisfied. Assuming 
that such variations exist, we establish the limit equation under the nondegeneracy condition that 
Q h approach Q along a direction U G Tq50(3) for which dur(Q) ^ 0. 
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Remark 1.3. Condition (jl.8j) of [16] is of technical importance. Notice that, in view of (11.20 
resulting in DW(F) = for all F G 50(3), ([L8]) is equivalent to: 

VF G M 3x3 < Cdist(F, 50(3)). 

Using the last assumption in (fl~2|) . the above implies that: \DW(F)\ < CW(F) 1 / 2 for all F G 
R 3x3 . Hence, roughly speaking, W has a quadratic growth and we see that ()1.8|) is actually 
quite restrictive. Independent from our research, Mora and Scardia [15] has presently established 
a result complementary to ours where the requirement (|1.8p is relaxed, while the equilibrium 
condition of (jl.3p is understood in a different manner, related to Ball's inner variations and the 
Cauchy stress balance law pp. 
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2. CONVERGENCE OF WEAK SOLUTIONS TO THE EULER-LAGRANGE EQUATIONS (EQUILIBRIA) 

OF THE 3D ENERGIES 

We first gather the relevant information from |10j : 

Lemma 2.1. [10] Let u h G W 1,2 (S h , M 3 ) be a sequence of deformations of shells S h . Assume 
ftl-4\) o-nd let the scaled energies E h (u h )/e h be uniformly bounded. Then there exists a sequence of 
matrix fields R h G W l,2 (S, M 3 ) with R h (x) G £0(3) for a.e. x G S, and there exists a sequence of 
matrices Q h G 50(3) such that: 

(i) \\(Q h ) T R h -U\\ m , Hs) <CVe^/h. 

(ii) h/\^e^((Q h ) T R h — Id) converges (up to a subsequence) to a skew- symmetric matrix field 
A, weakly in W 1,2 (S). 

Moreover, there exists a sequence c h G M 3 such that for the normalized rescaled deformations: 

y h (x + tn) = (Q h ) T u h (x + h/h tn) - c h 

defined on the common domain S h ° , the following holds. 

(iii) y h converge in W 1,2 (S h °) to it. 

(iv) The scaled average displacements V , defined in M.12\) converge (up to a subsequence) in 
W 1,2 (S) to some V £ V, whose gradient is given by A, as in $1.5]) . 

(v) h/Ve^ sym W h converge (up to a subsequence) in L 2 (S) to some B tan G B. 

The statements in Theorem 11.11 (i), (ii), (iii) are contained in the Lemma above. It therefore 
suffices to use the extra assumptions (ll.lOp and f 1 1 . 8 [> to recover equations (I1.13|) and (I1.14|) as 
h^O. 

We start by rewriting the equilibrium equation (|1.10|) in a more convenient form. Clearly, every 
variation (p h G W 1,2 (S h ,M?) can be by a change of variables expressed as: 

(2.1) 4> h {x + tn) =^(x + th /hn), 

for the corresponding ip G W 1,2 (S h °, M 3 ). Then, (jl.lOp becomes: 

n rho/2 

h 2 Ve h / f(x) f ip(x + tn) dt dx 

JS J-ho/2 



-ho/2 
ho/2 



(2.2, 

= h 

IS J-ho/2 



h f det(Id + th/h U)DW(Vu h (x + th/h n)) : V(f) h (x + th/h n) dt dx. 

JS J-ho/2 
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Notice also that: 

(2.3) V4> h {x + th/h n) = Vip(x + tn) ■ P(x + tn), 

where the matrix field P G L°°(S °, M 3 ) has the following non-zero entries: 

P(x + tn) tan = (Id + t/i//i n(x)) _1 (Id + tU(x)), n T P(x + tn)n = h /h. 
In view of Lemma EU define the matrix fields E h ,G h G L 2 (S ,/l °, M 3x3 ): 

£ h = -LzWQd + V^G h ), G h (x + tn) = -^= ({R h ) T Vu h (x + th/h n) - Id) . 



V e" V e" x 7 

With this notation, recalling the frame invariance of W in (jl.2p we get, for every F G R 3x3 : 



1 



:DW(Vu h (x + th/h n)) : F 



1 



=LW(i? ft (Id + Ve^G" 1 )) : F 

=LW(Id + \f^G h ) : (R h )F = R h E h : F. 



In particular, 
by V~ 
(2.4) 



becomes, after exchanging ip to (Q 1 ) \ip, using (|2.3p and dividing both sides 



rho/2 



ho/2 



Q h ip dtdx 



det(Id + t/i//i n) {Q h ) T R h (x)E n (x + tn) : V0 h (x + th/forO dtdz 



L™1 
II 

\ f ' det(Id + th/holl) (Q h ) T R h E h : [(V tan ip){ld + th/hoU)- 1 (Id + tn)] dtdx 
JsJ-h /2 L -Its 



feo/2 
ho/2 



/■ /-ho/2 

+ ho f det(Id + t/i//i n) ((Q h ) T R h E h n) d H ip(x + tn) dtdx, 

JS J -ho/2 

where Vj an denotes gradient in the tangent directions of T X S. The subscript TS stands for taking 
the 3x2 minor of the matrix under consideration, for example: 'Vtan^P = [V^]t5- 

Lemma 2.2. The sequence G h converges (up to a subsequence), weakly in 1? (S h ° ,M 3x3 ) to an 
L 2 (S h °) matrix field G, whose tangential minor has the form: 



(2.5) 



G(x + tn) 



tan 



K 

2' 



Btan ~ —(A 2 ) tan 



+ ± (V(Afi) - AU) tan . 
ho 



Moreover, if A1.8\) holds, then: 

(i) E h converges (up to a subsequence) weakly in L 2 (5 fi, °,R 3x3 ) to the matrix field E = C%G. 

(ii) The sequence (Q h ) T R h (x)E h (x + tn) converges (up to a subsequence) to E, weakly in 
L 2 {S h °,R ix3 ). 

Proof. The convergence of G h and the formula (|2.5p follow from Lemma 3.6 and Lemma 4.1 in 
|10j . Convergence in (i) is a consequence of Proposition 2.3 in [16], where the crucial role was 
played by the following equivalent form of the assumption (|1 .8j) : 

VFGM 3x3 \DW(ld + F)\ <C\F\. 

Finally, (iii) is an immediate consequence of (ii) in view of Lemma 12. II (i) and the boundedness of 
[Q h ) T R h in L^iS 110 ). ■ 
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Lemma 2.3. The matrix field E G L 2 (S h °, R 3x3 ), defined in LemmalKE (i) satisfies the following 
properties, a.e. in S h ° : 

(i) En = 0. 

(ii) E T = E, that is: E is symmetric. 

(iii) E tan (x + tn) = C 2 {x, G tan (x + tn)). 

Proof. To prove (i), one needs to pass h — > in (|2.4p and use Lemma 12.21 (ii) to obtain: 

(2.6) / -f 1 (e(x + tn)n) d n ip(x + tn) dtdx = 0. 

J S J -ho/2 ^ ' 

Now, any vector field (j) G L 2 (S h °, M 3 ) has the form <fi = d^tp, where ip(x+tn) = fth Q /2 ^{x+sn) ds. 
Therefore (i) follows from (12.60 . 

By frame indifference (|1.2f) and the fact that W is minimized at Id, it follows that DW{F) = 
for all F £ 50(3). It implies that for all H £ so(3) there holds C 3 H = 0, and so: E : H = C 3 G : 
H = C%H : G = 0, proving (ii). 

The assertion (iii) follows from E = C3G and the reasoning exactly as in the proof of Proposition 
3.2 pj. ■ 

A more precise information, with respect to that in Lemma 12.31 (ii) is given by: 

Lemma 2.4. There holds: 

(i) || skew E h \\ L1{sh0) < C\/?. 

(ii) lim — 1| skew E h \\ LP ^ S h ^ = 0, for some exponent p G (1, 2). 

Proof. By frame indifference flT2|) one has: = DW(F) : HF = DW(F)F T : H, for all F G M 3x3 
and all H G so(3) (since HF is a tangent vector to SO(3)F at F). We further obtain that 
DW(F)F T is a symmetric matrix. Apply this statement pointwise to the matrix field F = 
Id + Ve T G h : 

= -J= ( DW (Id + \/^G h ) (Id + Ve^(G h ) T ) - (Id + Ve^G h ) DW T (ld + v^O^ 
Ve' 1 v 

= E h - {E h ) T + v 7 ^ (E h {G h ) T - G h {E h ) T " 

Hence the claim in (i) is proved, as by Lemma l2.2t 

||sym (E h (G h ) T )\\ L i( S h ) < O||^ ,l || L 2( 5 h ) ||G fe || L 2( 5 h ) < C. 

Now, (ii) follows from (i) in view of the boundedness of E h in L 2 (S ho ), (jl.4p . and through an 
interpolation inequality: 

i|| skew E h \\ LP{sh0) < h\ skew E h \\ e L1 \\ skew E^ 6 < Cjh^ = C (Ve^/h 2 ^ 9 h 2e -\ 

where 1/p = 9 + (1 — 0)/2 and 9 G (0, 1). Clearly, the above converges to 0, when 9 > 1/2. ■ 

Introduce now the two matrix fields E, E G L 2 (S, M 3 ) given by the 0th and 1st moments of E: 

rho/2 rh /2 
E{x) = f E(x + tn) dt, E(x) = f tE(x + tn) dt. 
J-ho/2 J-ho/2 

It easily follows by Lemma 12.31 (iii) and Lemma 12.21 that: 

/ho/2 . k \ 

C 2 (x,G t an(x + tn)) dt = C 2 [x,B tan - -{A 2 ) tan ) , 
-h /2 V If 
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(2.8) 



Etn. n (x) 



rho/2 

i hi 

J-ho/2 



■>■• tG tan (x + tn)) dt = j±£ 2 (x, (V(An) - An) 



tan ) 



We will now use the fundamental balance (|2.4p and the above formulas to recover the Euler- 
Lagrange equations (|1.13p . (|1.14p in the limit as h — > 0. 

Proof of the first Euler-Lagrange equation (11.131) . 

Use the variation of the form: ifi(x + tn) = (j)(x) in f|2.4|> . divide both sides by h and pass to the 
limit to obtain: 

c r ho/2 r i 

0=lim/f det(Id + t/i//i n) (Q h ) T R h E h : [V tan ^(s)(Id + th/hoU)- 1 ] dtdx 

h ~*° J S J-ho/2 L ITS 

(2.9) = f 1° Ets . V ian ^(x) dtdx = f £ 2 (x,B tan - ^(A 2 ) tan ) : [V0(x)] tan dx 

JS J-ho/2 JS V I J 

C-2 (x,B tan - ^(A 2 ) tan ^ : sym V<f>(x) dx 

where we have used Lemma 12.21 (i), Lemma 12.31 and (|2.T|) . Therefore, by density of {sym V^} in 
the space B, (|1.13p follows immediately. 

Proof of the second Euler-Lagrange equation (j!.14|) . 

Let V G V and denote by A the skew-symmetric matrix field representing VV, as in (jl.5p . 

1. We now apply (|2.4|) to a variation of the form: ip(x + tn) = tAn(x). For simplicity, write 
n = Aft G W 1,2 (S, M 3 ). Upon dividing (|2.4p by /i and passing to the limit, we obtain: 

(2.10) 

rho/2 

det(ld + th/h U) (Q h ) T R h tE h : [V tan n(x){ld + th/hoYl)' 1 ] dtdx 



= lim 

h^O 



S J-h /2 



TS 



+ -r- I T {{Q h ) T R h E h n) r](x) dtdx 
h JS J-h /2 

r rho/2 

+ f (t tiacen + t 2 h/h detU)((Q h ) T R h E h n) n(x) dtdx 

Js J-ho/2 



-ho/2 

where we used the identity: 

det(Id + th/h U) = 1 + th/h trace U + t 2 h 2 /hl det II. 
The first term in (|2,10p , in view of Lemma 12.21 (ii) , Lemma 12.31 and (|2.8p , converges to: 



/ f tE TS : VtanV(x) dtdx = / E tan ■ (Vr)(x))t an dtdx 

JS J-ho/2 Js 

[ C 2 (x, (V(An) - AU) tan ) : (Vn(x)) tan dx. 
Js 



ho 
12 



In turn, the third term in (|2.10j) converges to 0. This is because (Q h ) T R h E h n converge weakly in 
L 2 (S ho ,R 3 ) to En = 0, by Lemma O (ii) and Lemma (i) . Summarizing, (|2.10p yields: 

(2.11) 

1 f rho/2 -i r 

lim - / f ((Q h ) T R h E h n) An dtdx = -— / C 2 (x, (V(An) - AU) tan ) : (V(An)) tan dx. 
h~>0 n Js J-ho/2 12 J s 
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2. Now, apply (j2.4j) to the variation ip(x + tn) = V(x), and pass to the limit after dividing 
both sides of 413} by h 2 : 

I f(x)-QV(x) dx = lim [ f(x) ■ Q h V(x) dx 
Js Js 

fho/2 



(2.12) 



lim 

h^O 



S J -ho/2 

h /2 -y 



\{{Q h ) T R h -ld)E h 
h 



TS 



+ 



LI 



S J-h /2 



TS 



A(x) TS (U + th/ho&dj n) 
A(x) T sO-d + th/ho&dj IT) dtdx 



dtdx 



:= lim [I h + II h ], 

where we used the definition of the adjoint matrix: 

det(Id + th/h U) (Id + th/hoU)' 1 = adj (Id + th/h U) = Id + th/h adj U. 
Notice that, by Lemma 12. II (ii) and (jl.4p . the matrix field: 

l/h((Q h ) T R h - Id) = (Ve h ~/h 2 )h/Ve h ~((Q h ) T R h - Id) 

converges to kA, weakly in W 1,2 (S) and hence strongly in L 2 (S). Hence, by the weak convergence 
of E h to E and the uniform convergence of (Id + th/ho&dj II) to Id, the first term of (|2.12p 
converges to: 

/ ° 7 (AE) TS :A ts = k[ (AE)ts :A TS = k[ (AE) : 
JSJ-h /2 JS JS 

-k [ E : (AA) = k f E tan : (AA) t 
Js Js 



lim Ih = k 

h-^0 



A 



(2.13) 



I tan 



Co I x, B, 



'tan 



'M 2 



I tan 



: {AA)tan dx, 



where we also have used Lemma 12.31 and (12 .7p . 

3. Towards finding the limit of IIh, in (|2.12|) . consider first the contribution of the tangential 
minors. By Lemma 12.41 (ii) and since A e LP(S ho ) for all p > 1, one observes that: 

-h /2 



i r c 

(2.14) lim - / j skew E? an : A^ 

h ~*° n Js J-i 



0. 



Hence: 



(2.15) 



lim / -f 

h^oJ s J_ 



h /2 1 



h /2 



_ rpn 
h tan " 



— lim 



-ho/2 



A(x)tan(ld + th/h adj n) 

ho/2 



ho h-*0j s J^ ho /2 



tE, 



tan 



■ — lim / 

ho h^o Js 



Etan '■ (Atan^l 



Han 



adj n 
1 

~12 



dtdx 



1 



— lim / Etan 

ho h^o J s 



A tan adj n 



C 2 (x, (V(An) - AU)tan) ■ (AU) tan dx, 



where we have used (|2.8p and Lemma 12.31 (ii) , combined with the following formula, which can be 
easily checked for A tan £ so(2): 

A tan adj n = -(AtanHf. 
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Further, by (|2TL|k 



lim 



II 



ho/2 1 



ho/2 



f ((£ h ) T ™)( W") didx 



-ho/2 



1 /* r 

lim - / / ((Q h ) T R h E h n)(Afi) 
h ^ n JsJ-h /2 



(2.16) 



+ lim 

ft 



Li 



ho/2 
h /2 

+ 2 lim 

h->0,/ g J-ft /2 



((Q h ) T R h -ld)(E h n 



h~\ 



(An) 



— ( skew £ )n 
h 



(An) 



^ / £ 2 (x, (V(An) - AU) tan ) : (V(In)) ton dx. 
J s 



Indeed, \/h((Q h ) T R h — Id) converges to weakly in L 4 (S) while ^4n E L A (S) and ft converges 
to weakly in L 2 (S). Therefore the second term in (|2.16p converges to 0. The last limiting term 
there vanishes as well, by Lemma 12.41 (ii) as in (|2. 14[) . 

Finally, we have: 

adjn)((i) T n) didx = 0, 



(2.17) 



lim — 

ft-»0 ho Js J -h /2 



n T tE h 



tan 



tan 



because (E h ) T n converges to weakly in L 2 (S) by Lemma 12.21 (i) and Lemma 12.31 
Adding now (|2TT5|) . (I2TT6|) and (I2TT71) we obtain: 

1 



(2.18) 



lim II h 

h-*0 



12 



C 2 (x, (V(An) - AU) tan ) : (V(An) - AU) tan dx. 



Together with (p7T2j) and (j2~T3l) . the formula (i2~18l) implies (fTTill . 



3. Convergence of critical points of the 3d energy functionals 

In this section we prove Theorem ll.2i Proceeding as in the proof of Theorem ll.il one needs to ex- 
change the expression J sh DW(Vu h ) : V(j) h by that of lim e ^ J s h \ [W(Vu h + eV(/) h ) - W(Vu h )] dz. 
As shown below, the error given by the difference of these two quantities, converges to as h — > 0, 
after an appropriate scaling by powers of h and Ve^ and along the variations <fi h used in the proof 
of Theorem ll.il 

We first prove a more general lemma, in which we derive the optimal asymptotic properties a 
sequence u h must satisfy in order that the conclusions of Theorem 1 1 . 1 1 hold true. These properties 
will later be established for the critical points (11.91) of the functional J h . 

Lemma 3.1. Assume \1.2\) and U.8\) . Let u h G W 1,2 (S h , M 3 ) be a sequence of deformations, 
satisfying: E h (u h ) < Ce h where the scaling e h is as in jl.4\ )- For every ip G W 1 ' 2 (S h °, M. 3 ), 
consider the rescaled variations <p h given by \2. 1\) and define the corresponding error terms: 

£ h {$) = [ DW(Vu h ) : V/ - / f h <f> h . 
Js h Js h 

Then all assertions of Theorem \1.1\ hold, provided that: 
(i) lim -—£ h (ifi) = 0, for all ip G W 1 ' 2 (S h °, R 3 ), 

h^0 \ ph 
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(ii) lim 



1 



-_£ h (ip) = 0, for all V of the form ip(x + tn) = <j){x), (j) G W ' (S,R 6 ), and all tp 



of the form ifj(x + tn) = tAn{x) with A given as in U.5\) for some V G V, 
(Hi) lim p=£h(?/>) = 0, for all ip of the form ip(x + tn) = V{x) with V G V. 



h-*0 h 2 Ve h 

Proof. The proof follows by a direct inspection of the proof of Theorem 11.11 Indeed, (i) is needed 
to derive ([231) . (ii) serves for getting (j2l)]) and (|2TL|) . through (pTTu]) . while (iii) implies ([2TT2]) . ■ 

Theorem 11.21 is now a consequence of the following observation: 

Lemma 3.2. Assume {LM> and (EH). Let u h G W 1,2 (S h ,M. 3 ) satisfy: E h (u h ) < Ce h with the 
scaling e h is as in |j.^[ ). // $1.9\) holds then the conditions (i), (ii), (iii) in Lemma \3.1\ are fulfilled. 

Proof. In view of (jl.9p . it is enough to prove that (i), (ii), and (iii) in Lemma [3.1 1 hold with £h{ip) 
replaced by a more convenient error quantity: 



4(V0 = iim n 4, e W, 

e— >U 



S h 



1 



ty(Vn h + eV0 h ) - W(ViA 



■LW(V« fe ) : V<£ fc dz, 



and the rescaled variations given by (|2.ip . 

Define the good sets: ft v = {z£ S*; dist(Vu h (z), 50(3)) < 5 and e|Vc^(z)| < 6}, with 5 > 
small enough for W to be C 2 in the open neighborhood {F G M 3x3 ; dist(F, 50(3)) < 35}. We will 
estimate £^ e (ijj) by writing it as a sum of two integrals: one over flh,e and the other over S h \fth,e- 
Apply the mean value theorem to the continuous function DW in the first integral, while in the 
second integral we use the assumption 



KM 



ds : V(f) h dz 



(3.1) 



[ [ \DW{Vu h + esV<p h ) - DW{Vu h ) 
Js h Jo 1 J 

<Ce I \Vcj) h \ dz 

+ C [ [ [dist(Vn h + esV(/) h , SO (3)) + dist(Vi/\ SO{3))} ds ■ |V^| dz 
Js h \n h e Jo L - 1 



< Ce 



[ \V(f> h \ + C [ dist{Vu h ,SO(3))\V(t> h \+Ce [ 
Ja hie Js k \n hie Js h \n K , 



h\2 



We see that the first and the third term above converge to as e — * 0. To treat the second term, 
notice that by (|1.2p . the energy bound, and (|2,3I) : 



(3.2) 



f dist(Vn /l ,50(3))|V^| <C f W(Vu h ) 1/2 \V(f> h \ 
Js h \n h ,, Js h \n h ,, 



< C 



1/2 



Observe also that: 



I |V/| 2 < C I 
Js h \n h « js 



/i|v tan ^| 2 + -|a^| 2 



dz, 
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where the set u>h,e = S h ° \{x + tn; x + th/hon G Its measure can be estimated as: 

K,| < c/h\s h \ n hje \ 

< C/hUz G S h ; dist(Vu h (z),S0(3)) > 5\ + \z G S h ; e\V<f>\z)\ > 5\\ 



(3.3) 



In particular: 



<C/h\ [ W{Vu h ) + e 2 ( \V<f> 
VJs h Js h 



(3.4) lim lim \ujh J = 0. 

We now prove (i). Passing to the limit in (|3.1|) and (|3.2|) . and using (|3.3|) with (|3.4p we obtain: 
lim lim (^)l < Climhm/i^HV^II^fcMo, i 



To prove (ii), consider: 

limlim^=|4 e (^)| < C limlimh- 1/2 \\V<t> h \\ L 2 (sh \ 



1 

Or— \ ' ™"' ' h 



< C lim lim ( IIV^n^H^^ ^ + rll%0ll 



The first limit above is by (|3.4p . Concerning the second term, it may be dropped for ip{x + tn) = 
4>{x), while in the other case when 4>{x + tn) = tAn(x) we have d^ip = An G W l ' 2 {S) and hence: 

lim limi||^||L2 Ke) < hm limCAI^I^II^ke^ < lim C/h{e h )^ = 0, 

in view of (jl.4p . 

To prove (iii) for ip(x + tn) = V(x), recall that VV G W 1 ' 2 (S,R 3 ) and write: 
lim lim — -j=\S' h M)\ < lim lim C/h\\VV\\ L 2 (u) , ) < lim lim C/h\w h J 1/3 || VV \\ L 6 (S) = 0, 
as before. This achieves the proof. ■ 

4. The limiting rotations Q 

In this section we will derive the third Euler-Lagrange equation (after the first two (|1.13j) and 
(11.14p ). corresponding to variation in Q G 50(3), and under certain nondegeneracy condition. 
We first notice that the limiting Q necessarily satisfies the constraint of the average torque: 

(4.1) r(Q) = I f x Qx dx = 0. 



Js 

The main difficulty arises now from the fact that the variations must be taken inside 50(3) in a 
way that this constraint remains satisfied. Assuming that such variations exist, we establish the 
limit equation under the additional condition that Q h approach Q along a direction U G TqSO(3) 
for which djjT{Q) / 0. 

In what follows, the crucial role is played by the function g(Q) = f s f ■ Qx dx defined on 50(3). 
Let B G K 3x3 be such that: g(Q) =B:Q, for all Q G 50(3). 
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Lemma 4.1. Assume the hypothesis of Theorem \l.l\ or Theorem Then the limit Q G SO(3) 
of Q h must satisfy: 



(4.2) / / • QFx dx = VF G so(3), 

Js 

or equivalently \J^. Another equivalent formulation of \J^.2^ is: skew(Q T B) = 0. 

Proof. First, for any given H G so(3), consider the variation <f) h = Hu h in the equilibrium equation 
(jl.lOp . Recalling that DW(Vu h )(Vu h ) T is symmetric (see the proof of Lemma \2Ah we obtain: 

(4.3) f f h -Hu h = [ DW(Vu h ) : HVu h = [ ( DW(Vu h )(Vu h ) T ) : H = 0. 
Js h Js h Js h v J 

Similarly, taking cf> h = ^(ex~p(eH)u h — u h ) in (|1.9|) . by frame indifference of W we get: 

f f h ■ Hu h = lim -I f h - {exp(eH)u h - u h ) = h lim - ( J h {e^{eH)u h ) - J h (u h )) = 0. 
J S h e^o e J S h c^o e V / 

Now, for any sequence of skew-symmetric matrices F h we have: 
j f-Q h F h V h = j J h ■ Q h F h {{Q h ) T u h -c h - id) 

(4.4) =^J sh fh - (Q hFh (Q h f) uh ~]^I sh fhdZ ' ~ ^ Jh ^ ' QhFhz dZ 

= I f-Q h F h x dx, 

Ve h Js 

where the first two terms in the second line above vanish by taking H = Q h F h (Q h ) T G so(3) in 
(|4.3p . and by the normalization of f h . Passing to the limit with h — > in (|4.4p . where F h = F, 
we see that: - J g f ■ QFV = lim h _> h/Ve 11 J s f ■ Q h Fx dx. This implies g2]). 

Clearly, (|4.2|) is also equivalent to = B : QF = Q T B : F for all F 6 so(3), which means 
exactly that Q T B is a symmetric matrix. 

To prove the other equivalent formulation of (|4.2p . notice that: 

I f ■ QFx = [ Q T f -Fx = -c F - [ Q T f xx = -c F [ f x Qx, 
Js Js Js Js 

where cp £ M 3 is such that Fx = cp x x for all i£K 3 . Since there is a one to one correspondence 
between vectors cp and skew matrices F, the proof is achieved. ■ 

Define now the set of the rotation equilibria: 

M = {Q G SO{3); skew (Q T B) = 0}. 

Our goal is to derive the third Euler-Lagrange equation, with respect to the variations of Q in 
M. For Q G M, let F £ so(3) be such that: 

QF = lim Qn ~ ® 



n-^oc ||Q n -Q||' 

for some Q n G converging to Q. Clearly, the above implies that: 
(4.5) skew (FQ T B) = 0. 
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Lemma 4.2. Under the hypothesis of Theorem \l.l\ or Theorem I LSI assume moreover that: 

Q h -Q 



lim 



t~o\\Q h -Q\\ 
Then for every F £ so(3) satisfying J^.5| ) there holds: 



QH, with skewiHQ 1 B) ^ 0. 



/ • QFV dx = 0. 



Proof. We will find a sequence F h S so(3), converging to F and such that, for all h: 



(4.6) 



/ • Q h F h x dx = 0. 



In view of (|4,4|) this will prove the lemma. Existence of such approximating sequence F h is 
guaranteed by the assumed nondegeneracy condition: skew (HQ T B) ^ 0. 
Firstly, notice that for Q h G M one can take F h = F. Otherwise, define: 

vh _ 1P (Q h fB : F 



skew {{Q h ) T B). 



|skew {{Q h ) T B)\ 2 



Then: 



/ / • Q h F h x dx = B: Q h F h = (Q h ) T B : F h 
Js 

(Q h ) T B : F 



and moreover: 
lim \F h -F\ 



lim 



= {Q h ) T B : F 
\(Q h ) T B:F\ 



skew ((Q h ) T B)\ 2 

\{Q h ) T B : F — Q T B : F 



(Q h ) T B : skew {{Q h ) T B) = 



lim 



h^o (skew {{Q h ) T B)\ h^o (skew ({Q h ) T B - Q T B)\ 



lim 



110* -Oil 



B : F 



/ 



skew 



Q h -Q 
\\Q h -Q\\ 



T 



B 



\H 1 Q 1 B :F\ 
(skew (H T Q T B)\ 



The last expression above equals to because of the nullity of its numerator: 

H T Q T B : F = Q T B : HF = Q T B : (HF) T = Q T B : FH = —FQ T B : H = 0, 
where we have used that Q T B is symmetric and (14. 5p . 
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